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Abstract 



t-h For the Oppenheim-Horodecki paradigm of nonclassical correlation, we propose a measure using truncations of a density 
^ matrix down to individual eigenspaces. It is computable within polynomial time in the dimension of the Hilbert space. 

^ 'keywords: nonclassical correlation, computational tractability 



\q 1. Introduction 

(N : 

Classical/nonclassical separation of correlations be- 
' c- 1 ' tween subsystems of a bipartite quantum system has been 
Q_ian essential and insightful subject in the quantum infor- 



^— > mation theory. The entanglement paradigm [ 1 , 2] is based 
^ on the state preparation stage: any quantum state that can- 
^ not be prepared by local operations and classical commu- 
, Q^ nications (LOCC) ^ is entangled. There are paradigms 
based on post-preparation stages, which use 
VP different definitions of classical and nonclassical correla- 
^ tions. On the basis of the Oppenheim-Horodecki defini- 
OO tion ASS], a quantum bipartite system consisting of 
^ subsystems A and B is (properly) classically correlated 
^ if and only if it is described by a density matrix having 
O a product eigenbasis (PE), p** = I*f k =\,\ e Jk\Vj)<yf\ ® 
^ jv^Xv^l, where d A (d B ) is the dimension of the Hilbert 
• • space of A (B), ejk is the eigenvalue of p A ® correspond- 
. ^ ing to an eigenvector |v^)® |v^). Thus, a quantum bipartite 
^ system consisting of subsystems A and B is nonclassically 
^ correlated if and only if it is described by a density matrix 
having no product eigenbasis. In this Letter, we follow 
this classical/nonclassical separation. 

This definition was introduced in the discussions by 
Oppenheim et al. |0, ^ on information that can be local- 
ized by applying closed LOCC (CLOCC, a branch fam- 
ily of LOCC) operations. The CLOCC protocol allows 
only local unitary operations and operations to send sub- 
systems through a complete dephasing channel. The clas- 



sical/nonclassical separation is linked to a localizable in- 
formation under the zero-way CLOCC protocol in which 
coherent terms are deleted completely by local players be- 
fore communicating under CLOCC. A bipartite state with 
a product eigenbasis carries information completely local- 
izable under zero- way CLOCC. The nonlocalizable infor- 
mation under zero-way CLOCC is a measure of nonclas- 
sical correlation. 

Other measures ifToL Hi . 12 , 13 ] than the nonlocaliz- 
able information were later proposed on the basis of the 
same definition of classical/nonclassical correlations. In 
particular, Piani et al. H 1 3TI designed a measure which 
vanishes if and only if a state has a product eigenbasis. 
It is in a similar form as quantum discord JH] and de- 
fined as a distance of two different quantum mutual in- 
formations that is minimized over local maps associated 
with local positive operator- valued measurements 11411 . A 
problem in a practical point of view is that the original 
nonlocalizable information and the Piani et al.'s measure 
both require expensive computational tasks to take mini- 
mums over all possible local operations. A similar diffi- 
culty exists in Groisman et a/.'s measure fioll for which a 
minimization is actually required to find a proper Schmidt 
basis (or dephasing basis) used to compute the measure in 
degenerate cases (namely, the cases where the eigenbases 
of one or both of the reduced density matrices of the state 
are not unique). A variant of Groisman et al. 's measure, 
the measurement-induced disturbance measure 1 15. fl^l . 
faces the same difficulty. 
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In our previous work [110 , an entropic measure G based 
on a sort of game to find the eigenvalues of a reduced den- 
sity matrix from the eigenvalues of an original density ma- 
trix was proposed. This measure can be computed within 



Preprint submitted to Elsevier 



March 29, 2010 



a finite time although it does not have a perfect detection 
range. Its computational cost is exponential in the dimen- 
sion of the Hilbert space. One way 11711 to achieve a poly- 
nomial cost is to introduce carefully-chosen maps similar 
to positive-but-not-completely-positive maps y, [l8|]. We 
pursue a different way in this letter. 

Here, we introduce a measure of nonclassical correla- 
tion for a bipartite state using the eigenvalues of reduced 
matrices obtained by tracing out a subspace after certain 
truncations of a density matrix. Its construction is rather 
simple as we see in Definition|2]of Section[3j The compu- 
tational cost is shown to be polynomial in the dimension 
of the Hilbert space. Although the measure is imperfect in 
the detection range and possesses no additivity property, 
it is practically useful as an economical measure invariant 
under local unitary operations. 

This letter is organized as follows. We begin with a 
brief overview of the measure G in Section |2] The mea- 
sure M is introduced and its properties are investigated in 
Section QJ An extension the detection range and the opera- 
tional meaning of M are discussed in Section^] Section[5] 
summarizes this work. 

2. Brief overview of the measure by partitioning eigen- 
values 

We first briefly overview the measure G, an existing 
measure computable in finite time. In the context of bi- 
partite splitting, it is defined as the minimized discrepancy 
between the set of the mimicked eigenvalues of a local 
system (say, subsystem A), {e;}- =1 , and the set of the gen- 

uine eigenvalues of the local system, {e,}- =1 . Here, e^s are 
calculated by (i) partitioning the d A xd B eigenvalues of the 
original bipartite state p AB into d A sets; and (ii) calculat- 
ing the sum of the elements for each set. The discrepancy 
in view from one side (from Alice's side in this context) 
is defined as 



F A (p AB ) = min 

parti tionings 



^(e, log 2 e t - a log 2 a) 



Similarly F (p ) is defined. The measure is defined as 

G(p AB ) = max[F A (p AB ) 5J F B (p AB )]. 

A drawback of the measure is that the number of com- 
binations of eigenvalues that should be tried in the mini- 

mizationis rf A rfB QBX (rf A_ 1)rfB QBX ••• x [t/ A^A„ 1)]rfB C> = 
(d A d B )l/(d B \) dA a 2 dAdBl0 ^ dA when the subsystem of 
concern is A [(d A d B )\/(d A \) dB a 2 rfArfBlo & rfB when it is 
B]. Indeed, this complexity is better in practice than that 



for minimization over all certain local operations required 
for calculating quantum deficit [6], quantum discord and 
Piani et. al. 's measure. The complexity of a minimiza- 
tion over all local operations for a subsystem, say A, is 
<9[poly(<A d B ) x 2 (rf ^ 2 '°g2 f ] with c the number of values 
tried for each parameter of a local operation. The com- 
plexity for computing G is smaller in the range d A ,d B < 
c. Computing G is, however, still very expensive. 

3. Measure based on partial traces of truncated den- 
sity matrices 

A measure that is computable within realistic time is 
desired for practical use. We introduce in the following 
a measure that achieves a realistic computational time, 
namely, polynomial time in the dimension of the Hilbert 
space. 

3.1. Introduction of the measure 

Let us begin with a basic definition. 

Definition 1. Let us write the eigenspace corresponding 
to the eigenvalue i] of the bipartite density matrix p AB as 
span{|v^)}^.j where dP is the dimension of the eigenspace 
and |v?)'s are the eigenvectors. Let us define a "truncated" 
density matrix down to the 77 eigenspace as 



k=i 



Proposition 1. Consider p I? introduced above. The eigen- 
values of the reduced matrix TrBp'' (Trpjfl ) of the system 
A(B) are integer multiples ofrj if p AB has a product eigen- 
basis. 

Proof. Suppose p^ has a product eigenbasis lla,)}^ X 
{\bj)}f =l . Then/5" becomes 77 2^ |a£><a>|^><^| where 
|fl^)'s (|^)'s) are some d v vectors, with possible multiplic- 
ities, found in {Ja,-)}^ ({1^)}^). Since |a,->'s (|fy>'s) are 
orthogonal to each other, it is now easy to find that the 
proposition holds. □ 

Remark 1. The sum of the eigenvalues of TrBp'' and that 
of TrA/5'' are both equal to Trp' ? = 77^. This property is 
tacitly used in the proofs we provide hereafter. 

Now we define the new measure of nonclassical cor- 
relation on the basis of Proposition Q] 
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Definition 2. For a given 77, consider the eigenvalues of 
the reduced matrix Tr^ff 1 . For each of the (dim Tr^p' 1 ) 
eigenvalues, the discrepancy from the nearest integer mul- 
tiples of 77 is calculated. Let us write the sum of the dis- 
crepancies as s 11 for a given 77. The measure of nonclas- 
sical correlation from the view of the subsystem A is de- 
fined as M A {p AB ) = Ya V J 1 - m tne same way, M B (p AB ) is 
defined. The new measure is defined by their average 

M(p AB ) = ^[M A (p AB ) + M B (p AB )]. 

The improvement in complexity is significant as we 
mentioned: this measure is computed by testing eigenval- 
ues of at most 2d A d B reduced matrices. The total com- 
plexity is dominated by the complexity of diagonalizing 
the original density matrix, 0(d A3 d B ), which is larger 
than the complexity of tracing out a subsystem for each 
truncated density matrix. In the exceptional cases where 
d A > d® 2 or d B > d Al , the complexity of diagonaliz- 
ing all the reduced matrices, 0[ma.x(d A4 d B ,d A d B4 )], be- 
comes the largest cost. Thus the total complexity, namely, 
the number of the floating-point operations in total to com- 
pute the measure, is 

0[max(J A V 3 , J A V, d A d BA )l 
Example. Consider the two-qubit state 



«■ = |(|oo><oo| + u+xi + 1) 



(1) 



with |+) = (|0> + |1»/ V2. This has no product eigenbasis 
while it is separable. The nonzero eigenvalue of g is 1 /2 
with the multiplicity two. The eigenvalue of TrB^ is 1/2 
with the multiplicity two; this leads to M A ($) = 0. The 
eigenvalues of Tr^g are (2 ± V2)/4; this leads to M B (g) - 
[2xl/2-(2+ V2)/4] + [(2- V2)/4-0xl/2] - 1- Vt/2. 
Therefore M(g) = (2 - V5)/4. 

Remark 2. For any pure state \<p) AB , M{\ip) AB {ip\) = if 
and only if \<p) AB is a product state. This is clear consider- 
ing the Schmidt decomposition of |^) AB . 

3.2. The measure is not perfect in the detection range 

We have achieved a reduction in the complexity of de- 
tecting a nonclassical correlation by introducing the mea- 
sure M. The measure does not have a perfect detection 
range as is expected from the fact that it does not test 
all the local bases unlike other expensive measures like 
(the minimized) quantum discord (with both-side test) or 



zero-way quantum deficit. For example, the measure M 
vanishes for the two-qubit state 

cr= l - (100X001 + 2|01><01| + 3|1+><1 + |) . (2) 
o 

Nevertheless, this state has no product eigenbasis because 
|0)(0| and l+X+l cannot be diagonalized in the same basis. 

Consequently, what one can claim is that a state for 
which the measure does not vanish is in the outside of 
the set B of the states for which the measure M vanishes, 
and hence in the outside of the set C of the states having 
a product eigenbasis, as illustrated in Fig. Q] Note that 
the set B includes some inseparable states while C does 
not. For example, the measure M vanishes for the state 
r, which is represented as a density matrix acting on the 
(3 x 3)-dimensional Hilbert space: 

r - \{\$)™(<P\ + |(A) AB <(AI + \0 AB (t\) 0) 



with 



\<t>) 

If) 



AB 



|0) A |1) B +|1) A |0) B 
V2 

AB _ |1> A |2) B +|2> A |1> B 
V2 

AB _ |2> A |0) B +|0> A |2) B 



V2 



This is because the nonzero eigenvalue of t is 1/3 with 
the multiplicity three and the eigenvalue of TrBT = TrAT 
is also 1/3 with the multiplicity three. The state t is in- 
separable because its partial transposition (/ ® T)t (here, 
T is the transposition map) has the eigenvalues -1/6, 1 /6 
and 1 /3 with multiplicities two, six and one, respectively. 
Thus it has been found that M cannot detect nonclassical 
correlation of this negative-partial-transpose (NPT) state. 
This example however does not weaken the measure M 
very much as there is a way to extend the detection range 
simply, which we will discuss later in Section |H 

3.3. Relative detection ability 

One can compare the measures M and G in their de- 
tection ability using the state a given in Eq. © and an- 
other particular state for a couple of qubits. 

The measure M is vanishing for cr while G is nonva- 
nishing. We have G{cr) = #(1/3) - H[{6 - Vl0)/12] a 
0.129 where H(x) = -^log 2 x — (1 - x) log 2 (l - x) is the 
binary entropy function. 

On the other hand, M is nonvanishing for the state 
cr' = 10X01/2 + (101X011 + I10X10D/4 with |0) = (|00) + 
|11))/ V2, while G is vanishing. We have M A (cr') - M B (cr') 
1/2 and M(cr') = 1/2. 

Therefore, M is neither stronger nor weaker than G in 
detecting nonclassical correlations. 
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Figure 1: The hierarchy of the quantum states. A: The convex set 
of the separable states. B: The nonconvex set of the states for which 
the measure M vanishes. C: the nonconvex set of the states having a 
product eigenbasis (C c A n B). Internal geometric structures of each 
set are not depicted, which should wait for future investigations. 



3.4. Invariance under local unitary operations 

The measure M is invariant under the local unitary op- 
erations ^^W^pw^® U B )p(U A ® U^. 

This is easily verified from the fact that (i) Tr^p 71 and 
Tr B CW A ® ( U B )(p'') have the same eigenvalues, and sim- 
ilarly, (ii) Tipfii and Tr A CZ/ A ® < i/ B )(p'?) have the same 
eigenvalues. 

3.5. Additivity properties do not hold 

The measure M has neither the subadditivity property 
nor the superadditivity property as we prove below. It is 
also shown to be not weakly additive. 

Here, let us denote a splitting of a system considered 
for the measure M by a subscript -|-. 

Proposition 2. Neither the subadditivity 

M A c|bd(p AB ® ^ CD ) < M A | B (p AB ) + M qD (^ CD ) 
nor the superadditivity 

M A c|bd(p AB ® o- CD ) > M A | B (p AB ) + M qD (^ CD ) 

holds. In addition, the weak additivity 
in 

M AlA2 ...| Bl B 2 ...(0p AiBi ) - mM Al | Bl (p AlBl ) 
i=i 

does not hold. 

Proof, (i) First we prove that subadditivity does not hold. 
Consider the state £ = cr^ ® cr CD with the state cr defined 
by Eq. ©. We have already found M A |b(o" AB ) = 0. Now 
we calculate M A c|bd(£)- The state £ has the eigenvalues 



e = 0, 1/36, 1/18, 1/12, 1/9, 1/6 and 1/4. Let us write the 
truncated density matrix down to the e-eigenspace as '. 
We have | 1/12 = (|001+>(001 + | + |1 + 00><1 + 00|)/12. 
This leads to Tr AC £ 1/12 - (|0+><0 + I + I + 0>(+0|)/12 
whose eigenvalues are 1 /24, 1 /8 and (with the multiplic- 
ity two). Similarly, £ 1/6 - (|011+><011 + | + |1 + 01)<1 + 
01|)/6 and Tr AC | 1/6 has the eigenvalues 1/12, 1/4 and 
(with the multiplicity two). For other £ e 's, Tr A c| e has the 
eigenvalues equal toe. Therefore, ^ B q BD (£) = (1 + 1+2+ 
2)/24 = 1/4. In addition, it is easy to find M A £ |BD (£) = 
because T^dI e has the eigenvalues equal to integer mul- 
tiples of e for every e. Consequently, M A c|bd(£) = 1/8, 
which is larger than M A | B (cr AB ) + M C |d(o~ CD ) = 0. This is 
a counterexample to subadditivity. 

(ii) Second we prove that superadditivity does not hold. 
Consider the state = ^ B ® g with the state g de- 
fined in Eq. £T|). We have already found M A |b(^" ab ) = 
(2- V2)/4. Now we calculate M A c|bd(£0- The state g has 
the eigenvalue 1 /4 with the multiplicity four and with 
the multiplicity 12. WehaveTr AC £' = (|0><0|+|+><+|)® 2 /4 
whose eigenvalues are (3 ±2 V2)/8 and 1/8 with the mul- 
tiplicity two. This results in M B ° BD (f ) = [3/4 - (3 + 

2 V2)/8] + (3 - 2 V2)/8 + 1/8 + 1/8 = 1- V2/2. In addi- 
tion, we have Tr BD £' = (|0><0| + |l>(l|)® 2 /4 whose eigen- 
values are 1/4 with the multiplicity four. This results in 
M ac|bd^') = °" Consequently, M AC |BD(f ) = (2- V2)/4, 
which is less than M m (g AB ) + M qD (f CD ) = 1 - V2/2. 
This is a counterexample to superadditivity. 

(iii) The above counterexamples shown in (i) and (ii) are 
also counterexamples to weak additivity. □ 

4. Discussions 

The main aim of introducing the measure M has been 
the computational tractability. Indeed, it is possible to de- 
cide whether a given density matrix has a product eigen- 
basis within polynomial time simply by diagonalization if 
there are only nondegenerate eigenvalues for the matrix, 
but it is, so far, not known to be possible if there are de- 
generate eigenvalues. Some existing measures, namely, 
Groisman et al. 's measure 1 10] and its variants, are com- 
putable within polynomial time if the dephasing basis is 
uniquely determined while it is not otherwise. The mea- 
sure M is, in contrast, always computable within polyno- 
mial time. 

As is expected for a tractable measure, M is imperfect 
in its detection range. The measure has been found to van- 
ish not only for the state © but for the NPT state ©. It is 
preferable for a measure of nonclassical correlation not 
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to vanish for obviously nonclassically correlated states. 
One may take the average of M and some entanglement 
measure to produce another measure of nonclassical cor- 
relation easily, which is nonvanishing for entangled states 
detectable by the entanglement measure. This approach 
was previously introduced in Ref. [ 17] to combine two 
different measures of nonclassical correlations. 

One can use any entanglement measure for this pur- 
pose. As an example, one may use the logarithmic nega- 
tivity lEl Up^) = log 2 [2Af(p AB ) + 1] with N(p AB ) the 
negativity II20L 12111 [namely, the absolute value of the sum 



of the negative eigenvalues of (/ ® T)p ]. The convex 
combination Mi = pM + (1 - p)L with < p < 1 is 
obviously stronger than M and L since they are neither 
stronger nor weaker to each other. It is, of course, invari- 
ant under local unitary operations. 

The evaluations of the measure M have been performed 
on its mathematical properties and computational cost. Let 
us discuss a rather conceptual problem. It is often of gen- 
eral interest to find an operational (i.e., process-based) 
meaning of a measure to justify the quantification. In this 
sense, the following argument may support the measure 
M. 

Suppose a dealer has a source generating quantum sys- 
tems that are in average described by the density matrix 
P = YjtjP' 1 - The source each time generates a system de- 
scribed by a state in an 77 eigenspace of p. The dealer each 
time sends Alice and Bob the corresponding parts of the 
generated system with the value of 77 attached. The dealer 
may modulate the source by a certain unitary operation 
common for all the times. The eigenvalues 77 of a density 
matrix p is unchanged under unitary operations. There- 
fore each p' 1 can be regarded as a data set carried in the 
"channel" 77, namely the eigenspace corresponding to 77. 

Consider which data can be secure data for Alice in 
the spectrum of Tr^p' 1 , against Bob's guesswork. The la- 
bel 77 of the channel is known. Therefore an eigenvalue 
written as an integer multiple of 77 is guessed by Bob within 
at most dP trials. In contrast, a guesswork to find an eigen- 
value that is not an integer multiple of 77 is computation- 
ally intractable for Bob: it takes exponential time in the 
size of the mantissa portion. The discrepancy between the 
eigenvalue and the nearest integer multiple of 77 is in other 
words the smallest shift from the insecure values. Recall- 
ing the definition, the measure M A can be regarded as the 
sum of the eigenvalue shifts for Alice over the channels 77, 
which are useful to protect her eigenvalues against Bob's 
guesswork. 

Similarly, M B can be regarded as the sum of the eigen- 
value shifts for Bob over the channels 77, which are useful 



to protect his eigenvalues against Alice's guesswork. The 
measure M has the meaning of the sum of the shifts aver- 
aged over the players. In this way, an operational meaning 
has been brought to M. 

The operational meaning introduced here is, however, 
presumed to be lost for the stronger measure Ml although 
it is controversial if the process of taking the average of 
different measures should be regarded as a nonoperational 
manipulation. 

5. Summary 

We have proposed an unconventional measure of non- 
classical correlation; its mathematical properties and its 
operational meaning have been investigated. It is invari- 
ant under local unitary operations while it is imperfect in 
the detection range and it has no additivity property. It is 
usable for a practical evaluation of quantum states because 
it is calculated within polynomial time in the dimension of 
a density matrix. 
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